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ABSTRACT 

m ■ 

Aims. Using a newly developed algorithm, we map, to the highest angular resolution allowed by the data, the temperature 
structure of the intra-cluster medium (ICM) within a nearly complete X-ray flux limited sample of galaxy clusters in 
the redshift range between z = 0.045 and z = 0.096. Our sample contains seven bright clusters of galaxies observed 
with XMM-Newton: Abell 399, Abell 401, Abell 478, Abell 1795, Abell 2029, Abell 2065, Abell 2256. 
Methods. We use a multi-scale spectral mapping algorithm especially designed to map spectroscopic observables from X- 
ray extended emission of the ICM. By means of a wavelet analysis, this algorithm couples spatially resolved spectroscopy 
with a structure detection approach. Derived from a former algorithm using Haar wavelets, our algorithm is now 
implemented with B-spline wavelets in order to perform a more regular analysis of the signal. Compared to other 
adaptive algorithms, our method has the advantage of analysing spatially the gas temperature structure itself, instead 
of being primarily driven by the geometry of gas brightness. 

Results. For the four clusters in our sample that are major mergers, we find a rather complex thermal structure 
with strong thermal variations consistent with their dynamics. For two of them, A2065 and A2256, we perform a 3-d 
^ ■ analysis of cold front-like features evidenced from the gas temperature and brightness maps. Furthermore, we detect a 

\Q 1 significant non-radial thermal structure outside the cool core region of the other 3 more "regular" clusters, with relative 

\^ \ amplitudes of about about 10 % and typical sizes ranging between 2 and 3 arcmin. We investigate possible implications 

of this thermal structure on the mass estimates, by extracting the surface brightness and temperature profiles from 
complementary sectors in the "regular" clusters A 1795 and A2029, corresponding to hottest and coldest regions in 
the maps. For A2029, the temperature and surface brightness gradients seem to compensate each other, leading to 
a consistent mass profile. For A 1795, however, the temperature structure leads to a significant mass discrepancy in 
the innermost cluster region. The third "regular" cluster, A478, is located in a particular sky region characterised by 
strong variations of neutral hydrogen column density, Nh, even on angular scales smaller than the cluster itself. For 
this cluster, we derive a spectroscopic Nh map and investigate the origin of Nh structure by discussing its correlation 
with galactic emission of dust in the infrared. 

Key words. Galaxies: clusters: general - Galaxies: intergalactic medium - X-rays: galaxies: clusters - Techniques: image 
processing - Techniques: spectroscopic 
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1. Introduction 



Forman fc Jonesl Il982l: ISlezak et all Il994l : iBuote fc Tsail 
1996t ISchuecker et al.ll200ll ). The current generation X-ray 



Clusters of galaxies are thought to form by accretion of less telesco (XMM-Newton, Chandra) now further allows us 

massive groups and clusters under the influence of grav- tQ accuratel and i nves ti g ate the temperature struc- 

ity. Following this scheme, they successively overcome some ture of thg JCM (gee | Ba uer et al.ll200l , provided that 

transient merging processes that alter their energy content ific al ithms hel at optimising a neC es S ary compro- 

and some relaxation phases, which lead to the almost viri- m . v ^ +1, ~e + ,, 

* ' mise between the precision 01 local temperature measure- 

ahsed systems we observe ... ments and the spatial resolution of temperature maps. 

the thermodynamical state 01 X-ray emitting mtra- 
cluster medium (ICM) depends on both the cluster merger For this reason, a number of spectral-mapping algo- 
history and some not yet perfectly understood physi- rlthms have been recently developed m X-ray astron- 
cal processes that drives its thermalisation, such as heat om y- All of them can be essentially subdivided m three 
conduction or viscosity. Revealed by X-ray observations, categories, depending on the general approach used. A 
the brightness and temperature structure of the ICM first approach combines X-ray imaging in several energy 
testifies to this final state. The brightness structure of bands > and colour analyses obtained from recombination 
the ICM has been extensively studied in the past (e.g. 01 multi-band information. The imaging part of these al- 
gorithms usually requires de-no ising by adaptive smooth- 
Send offprint requests to: H. Bourdin ing (e.g. Mark evitch et al.1 [20001. A second approach pro- 



2 



H. Bourdin and P. Mazzotta: ICM temperature structure within seven nearby and bright clusters of galaxies 



poses performing spatially resolved spectroscopy with re- 
quired signal-to-noise. To do so, the field of view is 
firstly sampled in independent spatial bins following a 
mesh refinement scheme, then spectroscopic estimations 
are performed within each bin. Several binning strate- 
gies have been proposed in t his context: adaptive binning 
(e.g. ISanders fc Fabianl [2QQlh , "conto ur binning" (ISandersI 
2006), or Voronoi tessellation ( e .g. ICappellari fc Copinl 
2003|). Proposed bv lBourdin et al.l (|2004f h a third approach 
consists in using Haar wavelet coefficients in order to couple 
a multi-scale spectroscopic analysis with a structure detec- 
tion scheme. 

Among these algorithms, only the last one performs a 
direct investigation of structure for the searched parameter 
itself -i.e. the ICM projected temperature- while the signal 
analysis is essentially driven by the geometry of gas bright- 
ness in all the other algor ithms. Here, we pres ent a new 
version of the algorithm of iBourdin et al.l ([2004), now im- 
plemented with B-spline wavelets and possibly enabling us 
to apply a more conservative thresholding strategy. Other 
improvements are related to adaptations to more recent cal- 
ibration of XMM- Newton instruments, which particularly 
allow us to combine data coming from the three European 
Photon Imaging Cameras (EPIC). 

This algorithm has been used to perform a systematic 
study of a nearly complete X-ray flux selected sample of 
seven clusters observed with XMM-Newton. Along this pa- 
per, we first discuss the sample selection and data prepara- 
tion in Sects. [2] and O then expose our data analysis scheme 
in Sect.JH and finally give a brief description of the thermal 
structure observed in each single cluster in Sect. [5l Before 
providing general discussions and conclusions in Sect. [UJ we 
discuss more specific issues related to thermal features ob- 
served within individual objects: isoradial thermal struc- 
ture in the relaxed clusters Abell 1795 and Abell 2029 
(Sect.[4j), cold fronts in the merging clusters Abell 2065 and 
Abell 2256 (Sect. ED . The additional Sect. El is related to the 
high neutral hydrogen column density variations observed 
across the field of view of Abell 478. 

The cluster radii are computed as angular diame- 
ter distances, assuming a A-CDM cosmology with Ho = 
70 km.s^.Mpc -1 , n = 0.3, n A = 0.7. 



2. The cluster sample 

St arting from th e X-ra y Brightest Cluster Sample (BCS) 
of lEbeling etTaD |l998), we selected a flux limited sam- 
ple of clusters. In order for the clusters to have an ob- 
served angular size close to the EPIC cameras field of 
view, we limited the cluster redshift selection to the range 
5z = [0.045, 0.096], corresponding to an angular size from 5 
to 10 arcmin. 

Our selection criteria returned eight bright clusters, 
namely Abell 2142, Abell 2029, Abell 478, Abell 1795, 
Abell 401, Abell 2256, Abell 399, Abell 2065, pub- 
lic XMM-Newton data being available for each of them. 
However, since the Abell 2142 observation is strongly al- 
tered by solar flares, we decided to remove it from our sam- 
ple and focus on the seven remaining clusters. 

Previous investigations, based primarily on the analysis 
of the X-ray surface brightness, classified 4 of the objects 
in our sample as major mergers and the remaining 3 as 
more "relaxed" clusters. The clusters are listed in Table [1] 



and are grouped according to this classification. The clus- 
ter redshifts and the ICM brightness from the BCS sur- 
vey are reported in the s ame table, with assoc iated P4/P0 
power ratio calculated bv lBuote fc Tsail ([1996) when avail- 
able. Expected to be an indicator of the dynamical evolu- 
tion stage of the cluster, this latter quantity is lower for 
the 3 most relaxed clusters than for the clusters classified 
as major mergers. 

3. Observations and data preparation 

3.1. Observations 

All data used for this investigation come from the EPIC 
XMM-Newton database. We use individual observations of 
Abell 399, Abell 401, Abell 478, Abell 1795, Abell 2029, 
and Abell 2065, and a multiple pointing observation of 
Abell 2256. For all of these clusters, data sets combine 
observations obtained from the three EPIC instruments: 
MOS1, MOS2 and PN. A summary of data sets is provided 
in Table [2j with associated exposure times. 

3.2. Data filtering 

In addition to the extended and optically thin emission 
of ICM, X-ray observations gather photon impact events 
related to other sources, such as spatially resolved X-ray 
emitting galaxies and the cosmic X-ray background (CXB). 
Moreover, observations may be transiently contaminated 
by solar flares. While the extended contributions of ICM 
and CXB are hardly separable, contributions from point- 
sources and solar flares can be isolated and removed from 
the observed signal through spatial and temporal wavelet 
analyses, respectively. 

In order to detect high solar flares periods and remove 
corresponding data sets, we analyse light curves with asso- 
ciated high energy events (10-12 k eV) and softer events (1- 5 
keV), respectively. As proposed bv lNevalainen et"aH (2005), 
this two-step analysis first enables us to isolate the most 
prominent flares at high energy, where ICM brightness is 
expected to be negligible, then detect some additional con- 
tribution of soft flares only. For each of these curves, we use 
a B3-spline "a-trous" wavelet algorithm in order to detect 
disruptions, and select the positive irregularities with am- 
plitude overcoming a 2<j significance threshold with regard 
to the light-curve fluctuation. This "cleaning" process has 
significantly lowered the effective exposure time of observa- 
tions, as reported for individual observations and pointings 
on Tabled 

In order to identify point-sources, we analyse EPIC- 
MOS event-lists which are more suited to imaging, and 
adopt an object separation algorithm der i ved fr om the 
multi-scale vision model of Biiaoui fc Rud (|!995f ). After 
detecting and imaging point-sources, we associate a mask 
to regions of the field of view dominated by point-sources 
contribution, and isolate events coming from these regions 
when analysing the signal. 

3.3. Data sampling - effective exposure 

XMM-Newton imaging cameras provide photon impact 
event lists with associated energy and position on the de- 
tector planes. In order to analyse our signal using discrete 
wavelet algorithms, we group events coming from various 
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Tabl e 1. Cluster sample with associate d redshift, coordinate s, total ICM flux from the BCS survey (Eb eling et al.l 
1998), and P4/P0 power ratio reported in Buote & Tsai (1996). Last column: average neutral hydrogen column density 
measured along line of sight of each cluster ([Dickey fc Lockmanlll990l ). 



Cluster 


Redshift 


Equatorial coordinates 


ICM flux 


Power ratio 


Neutral hydrogen column 


name 




(J2000) 


(1(T 15 W.m~ 2 / 10" 12 erg.s^.cm- 2 ) 


P4/P0 


density (10 24 m~ 2 ) 


Relaxed clusters 


A478 


0.088 


04 13 25.0 +10 27 54.0 


39.9 


0.025 


15.1 


A1795 


0.062 


13 48 53.0 +26 35 32.0 


68.1 


0.004 


1.2 


A2029 


0.077 


15 10 56.0 +05 44 42.0 


61.3 


0.050 


3.14 


Major merger clusters 


A399 


0.072 


02 57 53.0 +13 02 00.0 


28.8 




10.9 


A401 


0.074 


02 58 58.0 +13 34 00.0 


42.8 


0.157 


10.5 


A2065 


0.072 


15 22 42.0 +27 43 00.0 


22.2 




2.95 


A2256 


0.058 


17 03 58.3 +78 38 31.0 


49.5 


0.395 


4.1 



observations, if any, and sample them spatially into sky co- 
ordinate grids, with given angular resolution, a$. We also 
sample events in energy so as to perform spectral analyses, 
which leads to 3-D event cubes associated with each EPIC 
instrument, and sampled in position (k,l) and energy (e). 

For imaging and spectral-mapping purposes, we asso- 
ciate to these event cubes a set of local "effective exposure" 
E(fc, Z,e) and "background" cubes, nb c k Fbck(&,^e), with 
similar position-energy sampling, and use them for mod- 
elling the bolometric or spectroscopic ICM radiation; we 
report details about the background modelling in appendix 
lAl Let us define here the "effective exposure" E(fc, Z,e) at 
pixel [k,l] as the linear combination of CCD exposure times 
tccD(k,l,p) related to individual observations p, with cor- 
rection for spatial variations of the mirror effective area or 
so-called "vignetting factor", Aa m i rror (fc,/,e), transmission 
by other focal instrument, -i.e. reflection grating spectrom- 
eter (RGS)- trRGs, and detector pixel area with correction 
for gaps and bad pixels, accD For K observations, 

we get: 



scheme. Then, we locally estimate the gas temperature T 
and its fluctuation <tt by fitting a spectral model to the 
data, within each meta-pixel [fc,Z,j] of the different grids. 
A set of temperature maps T(k,l,j) is obtained, with as- 
sociated noise expectation maps <j^{k,l,j). Filtering the 
T(/c,Z,j) and <jt(^,Z,j) maps using high-pass analysis fil- 
ter enables us to code the temperature variations as wavelet 
coefficients Wx(fc, l,j) with expected noise aw T (^ an d 
to detect significant temperature structures as wavelet coef- 
ficients with amplitude overcoming a significance threshold 
depending on crw T (^,Z,j). Finally, we map the gas tem- 
perature using a Tikhonov regularised thresholding of the 
wavelet transform. 

More detailed in iBourdin et all ([2004), this general ap- 
proach has been adapted to the present study. In particular, 
the local spectral fitting is now implemented using an up- 
dated plasma emission code and now allows a multiple pa- 
rameter estimation. Furthermore, a B-spline wavelet trans- 
form is now used instead of the Haar wavelet transform. 
Below we describe both of these improvements in detail. 



K 

E(k,l,e) = y^ccp(M,p) x Aa mirror (^,Z,e,p) 
P =i 

x tr RGS (k,l,e,p) 

x accD(M,p) (1) 

4. Data analysis 

4.1. ICM temperature mapping 

In order to map the ICM temperature structure, we have 
built a spectral mapping algorithm coupling a spectroscopic 
and multi-scale analysis of the X-ray signal, to a wavelet 

mapping of the searched parameter s tructure. 

Mainly following the scheme of the lBourdin et all ([2004) 
algorithm -hereafter the 2004 algorithm- we first sample 
the field of view using redundant and square grids with 
typical size sj = 2-^ao, {j G [0, j max ]}, according to a dyadic 

1 Information about these instrumental effects are 
provided in the following XMM-Newton Current 
Calibration Files (CCF), corresponding to each obser- 
vation epoch: RGS*_QUANTUMEF, XRT*_XAREAEF, 
EMOS*_LINCOORD, EPN_LINCOORD, EMOS*_BADPIX, 
EPN.BADPIX 



4.1.1. Spectral fitting of local ICM temperature 

The local estimation of ICM temperature is performed by 
fitting a normalised spectral model, F evt (T, Z, Nh, e), to the 
data set associated with meta-pixel [k,l,j]. Combining con- 
tributions of ICM itself, nicM^O Ficm(T, Z, Nh, e), and 
"overall background", nbck Fb c k(^,Z,e), this model is sam- 
pled in photon energies, e, and depends on the ICM tem- 
perature T, metal abundances, Z, and neutral hydrogen 
density column along the line of sight, Nh: 

n evt (M) F evt (T,Z,Nh,e) 
= E(fc, Z, e) x nic M (fc, Ficm(T, Z, Nh, e) 

+n bck F bc k(M,e), ( 2 ) 

The unnormalised ICM contribution Ficm(T, Z, Nh, e) 
is modelled from plasma radiation flux 0icm(T, Z, e), fol- 
lowing the A strophysical Plasma Emission Code (APEC, 
ISmith et al.ll2Q0lf ) as a function of the gas temperature T 
and heavy elements abundance s Z, being set according t o 
the solar element composition of lGrevesse fc Sauvall (|1998| ). 
It is obtained by red shifting and distorting 0(T, Z, e) in or- 
der to take into account X-ray absorption by the galactic 
neutral hydrogen ft(Nh, e) with given column density along 
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Table 2. Effective exposure time of each EPIC XMM-Newton observation. In brackets: fraction of the useful exposure 
time after solar-flare "cleaning". 



Cluster 


MOS1 effective 


MOS2 effective 


PN effective 


name 


exposure time (ks) 


exposure time (ks) 


exposure time (ks) 


A399 


8.3 (58.6 %) 


4.6 (32.8 %) 


4.2 (45.2 %) 


A401 


12.0 (92.5 %) 


10.0 (77.4 %) 


2.8 (35.1 %) 


A478 


30.3 (44.2 %) 


38.0 (30.5 %) 


47.0 (43.6 %) 


A1795 


22.5 (45.8 %) 


23.5 (47.8 %) 


19.3 (45.5 %) 


A2029 


6.5 (37.5 %) 


7.7 (44.4 %) 


6.4 (50.6 %) 


A2065 


20.4 (61.0 %) 


19.9 (60.3 %) 


12.2 (55.8 %) 


A2256 (1) 


8.7 (62.1 %) 


8.5 (60.3 %) 


4.8 (54.8 %) 


A2256 (2) 


7.9 (48.5 %) 


8.6 (53.0 %) 


5.9 (50.9 %) 


A2256 (3) 


8.6 (47.3 %) 


8.9 (48.6 %) 


5.3 (26.3 %) 


A2256 (4) 


7.8 (47.8 %) 


7.1 (43.3 %) 


6.5 (58.0 %) 



Table 3. The quadratic B-spline analysis filters. 
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Fig. 2. The quadratic B-spline scaling function <p(x) and 
dual wavelet ip(x) in ID. 



the line of sight, Nh, following absorpti on parameters of 
Bal ucinska- Church fc McCammor] ([l992). Introducing the 
instrument area A(e), and convolving ICM radiation flux 
by the detector response in energy i?(ej3, we get the radia- 
tion power measured by the instrument for a given emitting 
source at redshift z: 



Fi CM (T,Z,Nh,e) 



A[e) x ft(Nh, e) x 



0icm [T,Z,(l + z)e] 



1 



(3) 



where instrument area A(e) takes into account the ef- 
fective area of mirrors drnirr or 

(e), the filter transmission 
tr futer (e) and detector quantum efficiency gcc£>(e)E|: 

A( e ) = QccD(e) x tr futer (e) X (^mirror (e). (4) 

In order to get a robust estimation of ICM tempera- 
tures, whatever the local statistics are, the spectral fit- 
ting is performed by maximising the log-likelihood function 
log/(T, Z, Nh, e) = logF evt (T, e^), where our spectral 
model is summed on all energy channels, e^. 

4.1.2. Wavelet mapping of ICM temperature structure 

In order to map the ICM temperature structure at scale j, 
we analyse the spatial correlation of temperature measure- 



2 Detector responses are tabulated within the XMM-Newton 
redistribution matrix files (RMF) 

3 Information about these instrumental effects are 
provided in the following XMM-Newton current cal- 
ibration files (CCF), corresponding to each observa- 
tion epoch: XRT*_XAREAEF, EMOS*_QUANTUMEF, 
EPN.QUANTUMEF, EMOS*_FILTERTRANSX, 
EPN.FILTERTRANSX 



ments using wavelet coefficients Wtj(&, I), computed from 
the temperature maps Tj(k,l). A trivial solution to this 
computation has been proposed in the 2004 algorithm: fil- 
tering the temperature maps Tj(k,l) using Haar high-pass 
analysis filters enables one to get a set of Haar wavelet co- 
efficients WT^jf^ 0- Indeed, the Haar wavelet is the dual 
function of the top-hat smoothing kernel Tlj(k,l), which 
may be applied to the searched map To(/c,/) for comput- 
ing the maps Tj (&,/), at scale j. However, thresholding 
Haar wavelet transforms usually generates square artifacts, 
in particular when analysing regular signals. 

Due to the expected smoothness of our signal, a better 
analysis of its spatial correlations ca n be provided by the 
more regular B-spline wavel ets; see Curry fc Schoenberd 
(|!947f ) for definition and e.g. iMallatl (|l998f ) for application 
to wavelet bases. Indeed, we do not expect any strong dis- 
continuities in our signal, due to joint effects of both the 
instrument PSF and the 2-D projection of the gas tem- 
perature structure. Since B-spline wavelets are the dual 
functions of the m degree B-spline interpolation functions 
obtained by (m+1) self-convolving of a top-hat smooth- 
ing kernel, they can also be used in our context; start- 
ing from the temperature maps Tj(k,l) and convolving 
them (m) times by a top-hat smoothing kernel Hj(k,l) 
provides a new set of smoothed maps Ts( m )j(k : Z), whose 
dual wavelet coefficients are m degree B-spline wavelet co- 
efficients WT,5(m),j(^5 0- Introducing the smoothing length 
a = 2 J at scale j, we get: 



-S(r, 



Uj(k,l) - l[ fe _a ?fc+ a_ 1 ] X l^_a^ + a_ 1 ], 

(M,j) = ILj(k,l)W*T(k,lJ). 



(5) 
(6) 



Like the Haar wavelet, the B-spline wavelets can be 
projected in bi-o rthogonal bases. W e ther efore adopted 
the shift- invariant Coif man fc Donohol (|l995h algorithm for 
computing the wavelet coefficients WT,s( m ),j(&, /), follow- 
ing a similar scheme as for the 2004 algorithm. To do so, 
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Abell 401 Abell 399 




-6 -4 -2 2 4 6 -6 -4 -2 2 4 6 

Arc minutes Arc minutes 

R.A.(O) 02 58 56.68 Dec(0) +13 34 36.0 R.A.(0) 02 57 49.26 Dec(O) +13 02 50.2 



Fig. 1. EPIC-XMM "soft" (0.5-2.5 keV) exposures images of the clusters in our sample. 



the smoothed maps T5( m )(fc, are convolved with a set 
of three high-pass analysis filters h(k)g(l), g(k)h(l) and 
g(k)g(l), associated with the B-spline wavelet, which leads 
to the three sets of wavelet coefficients Wx^^jj^ffc, I), 

W T) 5(m) )J >(M) and W T ,S(m),j,d(M) : 
W T ,S(m),i,/i(M) = M%(0* T S(m)(M, j), 

W T> s(m),j, v (M) = g(k)h(l)*T s{rn) (k,lJ), 
W T ,s(m),i,d(M) = g(k)g(l)*Ts(m)(k,lJ). (7) 

In order to re-construct the signal To(fc,Z), the wavelet 
coefficients are first convolved using the high-pass inverse 
filters, h(k)g(l), g(k)h(l) and g(k)g(l) associated with the 
bi-orthogonal wavelet analysis and added to each other, 
which allows us to average the redundant information. 
Thus, the risk occurring when re-constructing the signal 
from the thresholded wavelet transform is significantly low- 
ered. 

Here we decided to work with the quadratic B-spline 
wavelet (m = 2, see Fig. [2]), which is regular enough for our 



purpose but more compact than B-spline wavelets of higher 
degree. The quadratic B-spline analysis filters are reported 
in Table O 



4.2. ICM brightness mapping 



In order to map the ICM brightness structure, we have built 
an imaging algorithm where the local brightness Lx(fc,Z) 
is first estimated from a brightness model associated with 
pixel [k,l], then iteratively de- noised in the scale-space by 
means of a discrete Haar wavelet analysis. 

We estimate the local ICM brightness Lx(fc,Z) for en- 
ergy band Ae, by correcting the number of events n evt (A:, I) 
detected at pixel (k,l), from the expected "background" 
contribution, nb c k(fc?0? and exposure weighted instrument 
area, a(fc, I): 
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Abell 401 




-6 -4 -2 2 4 6 

R.A.(O) 02 58 56.68 Dec(0) +13 34 36.0 



Abell 399 




-6 -4 -2 2 4 6 

R.A.(0) 02 57 49.26 Dec(0) +13 02 50.2 



Fig. 3. ICM temperature maps of the clusters in our sample overlaid to the soft (0.5-2.5 keV) X-ray brightness isocontours. 
The isocontour levels are logarithmically equispaced by a factor \/2. 



a(fc, I) = f E(fc, Z, e) F evt (T , Z , Nh , e) de, (8) 

JAe 

n ev t(M) - n bc k(fc,0 



a(fc, Z) 



(9) 



where the weighted instrument area a(fc, Z) is com- 
puted for a fixed ICM emission model, F evt (T , Z , Nh , e). 
"Corrected" photon maps obtained from local estimator 
Lx(fc,Z) and for energy band 0.5 — 2.5 keV are shown in 
Fig.ffl 

The signal de-noising is then performed by thresholding 
Haar wavelet coefficients Wq [fc, Z] , associated with a multi- 
scale analysis of Lx(fc, Z) using redundant grids. Driven by a 
significance criterion of detected structure, this threshold- 
ing is set after estimating the probability density function 
(PDF) of noise wavelet coefficients, P(Wq[/c, I]). To do so, 
we deduce P(Wg[fc, Z]) from the PDF of noise wavelet co- 
efficients associated with local event counts, n evt : 



P(w L ~M) 



P(w nevt M) _ 

a[fc, Z] 



(10) 



which enables us to use an analytical fo rm of 

P(W nevt [MD, first proposed bv iBiiaoui fc Jammall (|200lh 
for Haar wavelet analyses of uniform Poisson noise. 
Introducing the modified Bessel function of order m, 

(n evt ), we get: 



P(W nevt )= E e- ne -/ m (n evt )(5(W nevt -m). (11) 



5. 2D temperature structure 

The ICM temperature maps of each cluster in our sample 
are presented in Fig. G2 overlaid to the relative 0.3 — 2.5keV 
brightness isocontours levels obtained by wavelet imaging 
(see details above). The emission model F evt (T , Z , Nh , e) 
assumed for computing brightness maps has been obtained 
by spectral-fitting of an overall cluster spectrum excluding 
the core region of relaxed clusters. 
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The temperature maps have been computed from 
spectral- fitting within the .7 — 7.5keV energy band with 
fixed values of the redshift, metal abundances, and neu- 
tral hydrogen column densities. The redshift values have 
been set to those in Table [TJ while the metal abun- 
dances Z have been deduced from the spectral fitting of 
F evt (T ,Z ,Nh ,e). For each cluster except Abell 478, the 
neutral hydrogen column densities, Nh , have been esti- 
mated by spectral fitting, similarly to the metal abun- 
dances. The values obtained are consis tent with mea- 
surements by iDickev fc Lockmanl ([l990f h in the 1 de- 
gree neighbourhood of each cluster centre (see Table [1]). 
For A478, the Nh has be en left as free parameter, since 
IDickev fc Lockmanl (|l990l ) measurements were inconsistent 
with X-ray spectroscopic estimations (see Sect. [8] for de- 
tails). 

The EPIC-XMM field of view has been sampled in 256 x 
256 pixels, and the wavelet analysis has been unfolded over 
six spatial resolutions, allowing detection of features with 
typical size ranging from about 12 arcsec to 3.5 arcmin. The 
thresholding of the wavelet coefficients has been performed 
following the lDonohol (11995) thresholding approach, leading 
to threshold levels that are always higher than the noise 
fluctuation a (typically from 1.75 to 2.5 a with increasing 
analysis scales, thresholds lower than 2a being related to 
the highest resolution details, namely the cool core regions 
of relaxed clusters only). 

In the following we give a short description of each clus- 
ter separately. For convenience we group the clusters fol- 
lowing the dynamical classification, used in previous work 
and based on the X-ray morphology, of relaxed and merging 
systems. 

5.1. Relaxed Clusters 
5.1.1. Abell 1795 

As observed in X-rays, Abell 1795 is known as an ellipti- 
cal and relaxed cluster, as predicted by the relatively low 
P4/P0 power ratios (see Tabled]). The XMM- Newton ob- 
servation of Abell 1795 confirms a globally elliptical sym- 
metry for this cluster. Nevertheless, the XMM- Newton X- 
ray image in Fig. [1] shows the presence of a sharp surface 
brightness variation, at about 1 arcmin to the south with 
respect to the cluster centre. From a Chandra observation, 
this sharp surface brightness feature has been identified as 
a cold front and interpreted as the result of the sloshing of 
core gas within the cluster potential well (Mar kevitch et aD 
This scenario is consistent with what can be observed 
on surface brightness contours in Fig. O revealing at the 
same time a strong compression of isophotal lines across 
the cold front, and a shift of the cluster brightness peak 
towards the south, with regard to large radii isophotes. 

The temperature map of Abell 1795 in Fig. [3] 
shows a global elliptical symmetry with a cool core. 
If we perform a radial analysis, it is consistent 
with temperature pr o files derived from XMM-Newton 



(Ta mura et aD 120011 lArnaud et all 1200 ll: llkebe et al l 
2004 1) and Chandra observations (lEttori et all |l002; 
Mar kevitch et al.ll200ll : IVikhlinin et all 12005). We notice, 



(r > 2 arcmin) the cluster also presents some significant 
non radial thermal structure. In particular we find that the 
gas temperature is colder in a small sector to the north 
(kT ~ 5keV) and hotter elsewhere (kT ~ 6keV). 

5.1.2. Abell 2029 

Similarly to Abell 1795, Abell 2029 was also known as a reg- 
ular and relaxed cluster (see also the P4/P0 power ratios 
in Table [1]). Except for its very cen tral region where X -ray 
filaments have been observed (e.g. IClarke et al.l 12004} ), no 
brightness nor temperature anisotropics have been previ- 
ously detected for this cl uster. 

As already noticed bv lVikhlinin et all (|2005l ), this clus- 
ter is projected near the galactic region of North Polar Spur, 
the X-ray spectral contribution of which has been modelled 
and added to the background model of Eq. (jA.ip . This con- 
tribution has been fitted simultaneously with cosmic back- 
ground, in the same external region of the field-of-view, 
which led to a two thermal component model (kTl=0.22, 
k T2=0.49, njgjnj2 = 2 .05), consistent with the analysis 



however, some anisotropics with regard to the overall ellip- 
tical symmetry of temperature structure. First of all, the 
cool core appears as being shifted towards the south, as 
coinciding with the cluster brightness peak. At larger radii 



of IVikhlinin et all (|200£ 

The XMM-Newton image in Fig. [1] confirms that 
Abell 2029 is globally regular and elliptically symmet- 
ric (see also Fig. [3|). The gas temperature map in Fig. [3] 
shows the presence of a cooler core (kT ~ 5.5keV) and 
a positive radial temperature gradient that extends up to 
1.5 arcmin from the cluster centre. If we perform a ra- 
dial analysis, the temperature map is c onsistent with pre- 
viously publish ed temperature profiles ([Lewis et al.l [2003: 
IVikhlinin et al.l 120051 ) , except in the very central region 
(r < 10 arcsec) , where the Chandra data analysis of 
iLewis et"aD (|2003l ) indicates lower temperatures decreasing 
to about 2 keV. 

Nevertheless, as observed for Abell 1795, even if the 
temperature structure is elliptically symmetric in the in- 
nermost regions, at larger radii the temperature map re- 
veals some non radial thermal features with typical tem- 
perature variations of c± IkeV. In particular, we detect a 
cold region at 5 arcmin to the southeast of the cluster cen- 
tre (kT ~ 6keV). 

5.1.3. Abell 478 

Known to be regular and relaxed, Abell 478 has been ob- 
served recently by the Chandra and XMM-Newt on tele- 
scopes ([Sun et al.l l2003l : iPointecouteau et al J 120041 ). While 
the Chandra observation has reveale d brightness sub struc- 
ture in the very centre of the cluster ([Sun et al.l l2003), both 
of these observations have shown a regular structure at 
larger radii. The elliptical symmetry and regularity of the 
surface brightness of Abell 478 is also evident from the X- 
ray cluster image in Fig. [1] and confirmed by the low value 
of the P4/P0 power ratios in Table [lj 

The gas temperature map of Abell 478 is shown in 
Fig. [3l Due to the strong hydrogen column density ap- 
pearing across the field of view (Nh = 15.1 x 10 24 m -2 ), 
the gas temperature map of this cluster has been com- 
puted from a spectral fitting process with a free Nh param- 
eter. Using this procedure the temperature map appears 



4 Normalisation discrepancies with Vikhli nin et al.l (|2005f ) are 
due to different modellings of the CXB, our model alrea dy in- 
cluding a thermal component at 0.204 keV (see Eq. (|A.1|) V 
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remarkably regular with elliptical symmetry. This map re- 
veals a cool core (kT c± 4keV) with 1 arcmin radius, a 
radial temperature increase up to a plateau at about 3.5 
arcmin from the cluster centre (kT ~ 7.5keV), then a de- 
crease up to 5 arcmin. This overall structure is in agreement 
with the shape of Cha ndra temperature profile published by 
IVikhlinin et al.l (|2QQ5l ). Nevertheless, the absolute temper- 
ature val ues from our temperat ure map are lower than the 
values of IVikhlinin et al~ (l2005f), and mor e consistent with 



the analysis of Pointecouteau et al. (|2004l ) ; using the same 



data a s here. IVikhlinin et al.l ( 20051 ) and ISanderson et al.l 



(2005) already noticed these discrepancies and argued that 
they are probably related to a combined effect of differences 
between instrument PSFs, and the complex temperature 
structure of the very central region in this cluster. 

From both the XMM- Newton brightness and temper- 
ature maps we may conclude that Abell 478 is the most 
regular of our sample. Nevertheless, we should note that 
this cluster is located in a particular region of the sky char- 
acterised by strong variations of neutral hydrogen column 
density, on angular scales smaller than the cluster size it- 
self. In this condition, the gas temperature variations are 
more difficult to detect for this cluster than for the other 
clusters of our sample (see also Sect. [8] below). 

5.2. Merging clusters 
5.2.1. Abell 399 - Abell 401 

Although, independently selected in our cluster sample, 
Abell 399 and Abell 401 form a close pair separated, in 
projected distance, by 3 Mpc (36 arcmin). The brightness 
contours of both Abell 399 and Abell 401 (Fig. [3| are 
elongated along a N-N-E / S-S-W axis, which is the major 
direction of the pair. The Abell 401 brightness contours ap- 
pear mildly disturbed with a centroid shift to the northeast. 
Consistent with the expectation of the P4/P0 power ratios 
reported in Table [IJ they are less regular than the contours 
of Abell 1795, Abell 2029 and Abell 478. Abell 399 looks 
even more irregular; in particular a sharp edge can be ob- 
served to the southeast of the cluster core. 

Using the same data set as used in this work, 
Sak elliou fc Ponmanl (|2004l ) studied the gas brightness and 
temperature in specific angular sectors of both clusters and 
along the major direction of the pair. From this analysis, 
they conclude that, currently, the clusters are just starting 
to mildly interact and that the sub-features found in their 
inner regions are related to the individual merging histories 
of each cluster separately, rather than to the remnant of a 
previous merger of the two systems. 

The temperature maps of Abell 399 and Abell 401 are 
reported in Fig. [3] to the same spatial scale. We notice a 
higher average temperature for Abell 401 (kT c± 8keV) than 
for Abell 399 (kT ~ 7keV) , consistent with the analysis of 
Sak elliou fc Ponmanl (|2004t ). The temperature maps of both 
clusters appear as strongly irregular, without any ellipti- 
cal symmetry or large scale correlation with gas brightness 
structure, contrary to what is observed for the relaxed clus- 
ters in our sample. Moreover, instead of cool cores, their 
central regions host some hot substructure embedded in 
the colder ICM. W e further confirm the tempe rature incre- 
ment described by Sak elliou fc Ponmanl (|2004t ) to the south 
of Abell 401, in the direction of Abell 399. However, we do 
not detect any temperature increment toward the north of 



Abell 399 along the expected interaction axis of the clus- 
ter pair. More generally speaking, most of the temperature 
irregularities detected here do not have any morphological 
link with the interaction axis of the cluster pair, and have 
typical sizes more related to the central region of each clus- 
ter than to the overall cluster system. 



5.2.2. Abell 2065 



Abe ll 2065 has been previously studied by Mar kevitch et al.l 
(11999ft using ROSAT and ASCA data, and by 
Chatzikos_et_aL| (2006) using Chandra data. These 
observations revealed the asymmetric morphology of this 
cluster characterised by a highly elongated inner region 
that seems to link together the two cD galaxies located 
in the cluster centre. The temperature structure is also 
highly asymmetric with a hotter region to the southeast. 
Furthermore, the Chandra observation revealed two cold 
cores coinciding with the cluster central cD galaxies. 

The XMM-Newton temperature map of Abell 2065 in 
Fig. [3] shows the presence of a hot (kT ~ 6.5keV) bow- 
like region, from about 3 to 6 arcmin to the southeast of a 
cold cluster core (kT ~ 4.5keV). This feature by itself is not 
isothermal but rather appears to embed an even hotter sub- 
feature (kT 8.5keV) to the eastern cluster outskirts. The 
map also shows a cold region (kT ~ 4keV) located at about 
7 arcmin to the northwest of the cluster core, appearing to 
be linked to this core by an extended tail. 



Al ready reported in pre vious works (M arkevitch et al.l 

1999; Chatzi kos etHll200d ). the hot bow-like region is lo- 
cated next to an abrupt variation of gas surface brightness 
visible in Fig. [TJ For this reason, it probably indicates the 
presence of a cold front, as discussed in Sect. 



5.2.3. Abell 2256 

The X-ray image of Abell 2256 in Fig. [1] shows a complex 
structure with two X-ray peaks, one of which is coincident 
with the cluster central dominant galaxy, while the other 
is located at 2 arcmin to the northwest. This cluster has 
been observed with RO SAT and ASCA (e.g. iBriel et all 
119911 : lMarkevitchlll996l ) and more recently with Chandra. 
This latter observation revealed an even more irregular X- 
ray morphology, the presence of a third subgroup to the 
east, and a sharp edge to the southeast of the northwest 
peak that has been identified with a "cold front" (|Sun et al.l 
I2002D . 

The XMM-Newton temperature map of Abell 2256 in 
Fig. [3] shows a bimodal temperature structure along the 
cluster major elongation axes. Consistent with Sun et al. 
(2002), we find that the gas in the northwest peak, is the 
coldest of the cluster with a temperature of kT c± 4.5keV. 
Unlike previous work, our temperature map also shows a 
clear, hot (kT ~ 9keV), bow-like region to the east. This 
hot region is located just outside an abrupt variation of the 
gas surface brightness, also clearly visible from the compres- 
sion of the isocontour levels in Fig. OS As shown in Sect. 
this new feature is likely to be related to another cold front 
in the cluster ICM. 
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Fig. 4. Spectroscopic temperature (top), brightness (middle) and derived mass profiles (bottom) corresponding to the 
selected red and blue sectors of Abell 1795 (left) and Abell 2029 (right), respectively. Dashed lines on temperature and 
brightness profiles correspond to fits of the projected functions {T x (r), £ x (r)} used for deriving masses, see Eq. (|T5|) 
and ([T5|) . Error bars on temperature profiles are 68 % confidence levels. Dispersions of brightness and mass profiles are 
distribution variances. 



6. Radial temperature structure and cluster mass 

X-ray emission from round relaxed clusters of galaxies is of- 
ten used to estimate the cluster mass, assuming hydrostatic 
equilibrium and spherica l symm etry. Using hydro N-body 
simulations, iRasia et al.l ((2006) investigated the accuracy 
of the mass estimate and found that, in the best of cases, 
there is at least a 10% discrepancy between the true and 



the estimated mass. They claimed that one of the reasons 
for such a discrepancy is related to small, non radial ther- 
mal substructure, which seems to be always present, even 
in the most relaxed systems produced by any hydro N-body 
simulation. 

To test whether this problem may also be present in 
real clusters, we used our temperature maps and selected 
specific cluster sectors in which the gas temperature is ei- 
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ther hotter or colder than the gas mean temperature. From 
these sectors, we extracted temperature and brightness pro- 
files, then used these profiles to estimate the relative cluster 
mass profiles. 

For this test we consider only two of the three more 
relaxed clusters in our cluster sample: Abell 2029 and 
Abell 1795. Abell 478, the third round cluster in our sam- 
ple, was not considered for such a test because, as explained 
in Sect. [HI it is located in a particular sky region charac- 
terised by strong angular variations of the neutral hydrogen 
column density, Nh, across the cluster field of view. If not 
treated properly, these strong Nh variations may introduce 
spurious thermal features that cannot be easily disentan- 
gled from the real ones. 

ICM surface brightness and temperature profiles 

Brightness and temperature profiles associated with hot 
and cold regions of temperature maps for Abell 2029 and 
Abell 1795 are shown in Fig. [H with different colours. 

The brightness profiles S x (r) have been obtained by av- 
eraging the brightness estimator of Eq. ([9]) within logarith- 
mically equispaced elliptical annuli of N pixels {fc, Z} r , with 
normalisation by the exposure weighted instrument area 
a(fc,/) of Eq. ©: 



{k,l} r 

N ^ 



{k,l} r 



nevt(fc,l) - n bc k(fc,0 
a(fc, I) 



(12) 



We then estimated the brightness variance by consider- 
ing uncertainties to be only related to the Poisson fluctua- 
tion of n evt (/c, I): 



\ 



£ 

{k,l} r 



<Jk,l) 
a 2 (k,l) 



\ 



E 



n e vt(fc,Q 
a 2 (k,l) ' 



(13) 



The spectroscopic temperature profiles T x (r) have been 
estimated by fitting the spectral model F ev t£ \ (T,Z,Nh,e) 
of Eq. ([2]) to the data associated with a set of logarith- 
mically equispaced sector annuli. Given the high statistics 
available for each annulus, we were able to perform spec- 
tral estimations following a y 2 minimisation process, and 
get straightforward estimations of confidence intervals. 

Cluster mass profiles 

The surface brightness and temperature profiles extracted 
from the hottest and coldest sectors of A1795 and A2029 
(see Fig. |4]) have been used to estimate the cluster mass 
profile. To do so, w e adopt the approach proposed by 
IVikhlinin etHI {2006), that consists in modelling the 3-d 
density and 3-d temperature profiles and fitting the pro- 
jected quantities to the corresponding data set. 

The 3-d gas density is modelled by a double and modi- 
fied /3-model with 9 parameters. Introducing n p and n e , the 
proton and electronic density, respectively, we get: 



[n p n e ](r) = n 2 , 



(r/r c )~ a 



[1 + (r/r c ) 2 p-«/ 2 [1 + (r/r s ) 3 ] e / 3 



L 02 



[l + (r/r c2 ) 2 p 2 - 



(14) 



The projected surface brightness profile, E(r), as ob- 
served within a given energy band AE (here 0.7-2.5 keV), is 
then obtained by integrating the ICM brightness, cicm(T), 
along the line of sight: 



Sx(r) = 



d 2 (l + z) 4 



eiCM[T(n)] [n p n e ](n) dZ, (15) 



where z is t he reds hift of the source located at distance 
d, where r\ = VV 2 + Z 2 , and where cicm(T) can be modelled 
from the source radiation power of Eq. (|3J), so as to account 
for instrument response and effective area. Normalising 
with the same average quantities, T , Z , Nh , as used for 
computing the weighted instrument area in Eq. (|5J), we get: 



cicm(T) 



/ A ^Fic M (T,Z,Nh,e) de 
J A£; FicM(T ,Z ,Nh ,e) de' 



(16) 



The 3-d temperature profile is modelled by a 5- 
parameter broken power-law, with transition region: 



T(r) 



(r/n) ° 



°[l + (r/r t ) b ] c / b ' 



(17) 



then integrated along the line of sight to get a projected 
profile of "spectrosco pic-like" temperatures , as defined in 
Eqs. (14) and (15) of iMazzotta et all (|2004f ): 



Tx(r) 



/ w(ri) dl 



J w(n)T(n) dl, (18) 



with weighting factor w(l) = 

After fitting the projected ICM 
temperature profile T s (r) 



brightness S x (r), 

and temperature profile T s (r) to the observable set 
{E x (r), T x (r)}, and estimating the related 3-d density p(r), 
and temperature T(r), we use hydrosta tic equil i brium of 
the ICM to derive mass profiles (see e.g. ISarazinl (|l988f ): 



M(r) 



-3.68 x 10 13 M o T ( r)r 



d log p d log T 



d log r d log r 



(19) 



Best fits of the projected ICM brightness and tempera- 
ture profiles with associated cluster mass profiles M(r), are 
shown in Fig. [4j The confidence intervals on mass profiles 
have been estimated by minimising the distance between 
the projected models and a set of random realisations of ob- 
served profiles. These profiles have been obtained assuming 
Gaussian statistics around observed values {£ x (r), T x (r)}, 
with the constraint of rejecting realisations leading to the 
non-physical solution of non-monotonically increasing mass 
profiles. 
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6.1. Abel 1 1795 

The Abell 1795 ICM temperature map reveals a cold re- 
gion to the north of the cluster, while the gas temperature 
is observed as hotter elsewhere in the 2 to 5.5 arcmin range 
of cluster radii. We selected two complementary "cold" and 
"hot" cluster sectors in order to compute surface brightness 
and temperature profiles in Fig. 2J As we want to ignore ef- 
fects of gas thermal variations near the cluster core, we ex- 
cluded this region from our analysis and computed profiles 
within a radii range of 1.5 — 8 arcmin. We further centred 
sectors so as to fit the cluster brightness isophotes at large 
radii, therefore the cool core appears as being shifted with 
regard to our sectors in Fig. [H 

Consistent with what is observed on the temperature 
map, we notice that the absolute temperatures in the two 
sectors are significantly different. Furthermore, we observe 
that the two temperature profiles show different shapes: 
the profile corresponding to the coldest sector peaks at ap- 
prox 4 arcmin, while the complementary hottest profile is 
much flatter. The two complementary brightness profiles 
also show different shapes. 

It is worth noticing that the resulting mass profiles show 
significantly different shapes (see Fig. [4|). In particular, as 
for the temperature profiles, the mass profile of the coldest 
sector shows stronger gradient variation than the hottest 
one. These different shapes result in the observation that 
while the relative mass estimates are consistent with each 
other at large radii (r > 250 kpc) -a nd also consi s tent w ith 
previously publ i shed profiles of e.g. Ilkebe et all (|2004f ) or 
IVikhlinin et al.l (|2006[ )- they become significantly different 
towards the innermost part of the cluster (r < 250 kpc). 

We should conclude that, within the radii range inves- 
tigated, the assumption of hydrostatic equilibrium for this 
cluster may be not valid. 



6.2. Abell 2029 

The Abell 2029 ICM temperature map reveals a non sym- 
metric temperature structure within a region ranging from 
2 to 5 arcmin cluster radii. We selected the coldest region 
to the southeast for computing a first set of brightness and 
temperature profiles, and a complementary hot sector for 
computing the second set. 

As with Abell 1795, the discrepancy observed between 
temperature values on profiles in Fig. [H is consistent with 
the radial thermal structure of the temperature map. In 
particular, the profile corresponding to the cold sector is 
almost fiat (kT ~ 6.5keV), while the profile corresponding 
to the hot sector has a positive gradient up to approx 4 ar- 
cmin radius (where it reaches kT ~ 8.5keV), then decreases 
down to kT ~ 6keV at larger radii. It is important to say, 
however, that due to the limited statistics available for fit- 
ting temperatures in the coldest region of A2029, the tem- 
perature discrepancies we find are only marginal for most 
of the radial bins. However, as with A1795, we find that 
the surface brightness profiles of the two complementary 
sectors have different shapes. 

Despite this, and in contrast with what we find for 
Abell 1795, the mass profiles derived for the two com- 
plementary cluster sector seem to have remarkably sim- 
ilar shapes that lead to a consistent mass estimate over 
the entire radial range considered (r G [200,800] kpc). 



This estima t e is a lso co nsistent with pr e vious works of e.g. 
iLewis et al.l (|2003f ) and IVikhlinin et all (|2006t ). 

For this cluster, we can say that the temperature and 
surface brightness gradients compensate each other very 
well, leading to the same mass estimates regardless of the 
sector used. Thus, in this specific case, we may conclude 
that the hydrostatic equilibrium assumption is well satis- 
fied. 



7. Cold fronts 

The temperature maps of both merging clusters Abell 2065 
and Abell 2256 show two hot bow-like regions to the south- 
east of each cluster centre, which seem to be located next 
to abrupt variations of the gas surface brightness, as visible 
on photon maps in Fig. [1] or isocontour levels in Fig.[3j For 
this reason, they are likely to be related to cold fronts sep- 
arating the dense and moving cluster cores from the hotter 
ICM of the cluster outskirts. 

In order to investigate this hypothesis, for clusters 
Abell 2065 and Abell 2256 we extracted the ICM surface 
brightness and temperature profiles corresponding to sec- 
tors shown in Fig. [5j Located across the brightness front 
regions, these sectors follow the brightness isophotes and 
match regions with almost uniform isoradial thermal struc- 
ture. 

As expected from temperature and brightness maps, we 
clearly see two temperature jumps on the profiles located at 
a jump radius r 3 ;, also corresponding to a change of slope in 
surface brightness profiles. In order to check whether these 
features can be related to jumps in the 3-d distributions of 
the gas density and temperature, we modelled these distri- 
butions by two disrupted functions. 

The gas density profile is modelled by two independent 
/3-models corresponding to regions located inside and out- 
side the jump radius r^, respectively: 



[n P n e ]cf(r) 



D?nl(r/rj) 2 



l+(rj/r cl f 
l+(r/r c i) 2 

l+(fi/r e2 ) 2 " 
l+(r/r c2 ) 2 



3/3i 



3A 



, Vr G [0, Tj[ 

, Vr G [r*j, oo [ 
(20) 



while the gas temperature is modelled by a step func- 
tion: 



T c /(r) 



T , 
D T T 



Vr G [0,rj[ 
Vr G [r j , oo [ 



(21) 



The 3-d distributions [n p n e ] c f(r) and T c /(r) are pro- 
jected according to Eq. (fT5|) and (fT8|) in order to derive 
X-ray brightness and "spectroscopic-like" temperature pro- 
files, S x (r) and T x (r), which enables us to estimate all free 
parameters of Eq. ([20]) and ([2T]h including r^, by fitting the 
projected models to the data. Best fits of X-ray brightness 
and "spectroscopic-like" temperature profiles are shown in 
Fig. O with associated 3-d distributions of gas density and 
temperature. Similarly to the derivation of mass profiles in 
Sect.[6j the confidence intervals on 3-d profiles have been es- 
timated by minimising the distance between the projected 
models and a set of random realisations of the observed 
profiles. 
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Fig. 5. From top to bottom: gas brightness, spectroscopic temperature, and derived density and temperature profiles 
corresponding to sectors shown on the top maps for Abell 2065 (left) and Abell 2256 (right), respectively. Dashed lines on 
temperature and brightness profiles corresponds to fits of the projected functions {T x (r), S x (r)} (see Sect.[7j)- The front 
modelling and fitting region is bounded by vertical dashed lines on the gas brightness and spectroscopic temperature 
profiles, while the fitted front position rj is reported by a vertical dot-dashed line. Error bars on the temperature profiles 
are 68 % confidence levels. Dispersions reported on the 3d profiles correspond to variances on each distribution. 
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7.1. Abel I 2065 

The Abell 2065 sector is located within 5 arcminutes to the 
southeast of the cool elongated central cluster region visible 
in Fig. [TJ 

Fitting the disrupted density and temperature profile of 
Eqs. ([201) and (EB) for the southeastern sector of Abell 2065 
leads to a jump radius of rj ~ O.mlo^ 

Mpc, a density 

jump factor between regions located immediately above and 
below r h of Dj ~ 1.29±Hi> 

and a temperature jump factor 

of Dt — 1.49+q 14. These values yield an almost continuous 
gas pressure across the front. Interestingly, we notice that 
the value of rj is consist ent with the loca t ion o f a density 
jump already reported by Chatzi kos et al.l (|2006l ), who per- 
formed a sector analysis of the same cluster region using 
Chandra data. The additional detection of a temperature 
jump using XMM- Newton enables us to identify the feature 
as a cold front. 

7.2. Abell 2256 

The profiles shown in Fig. [5] correspond to a sector located 
within 8 arcminutes to the southeast of the eastern cluster 
peak. As with Abell 2065, we detect a cold front feature on 
the profiles. The front is located at Tj ~ 0.240±q.oo2 M P C 
from the sector centre, with density and temperature jump 
factor of Dj = 1.30i™2 and D ? = ^ -o.it respectively. 
Also here, the values of jump factors are consistent with 
continuous gas pressure across the front. 

7.3. Possible origin of the fronts 

As shown by IChatzikos et al.l (|2006l ) using Chandra data, 
the cold front feature seen in Abell 2065 is located to the 
southeast of a dual cluster core, which is likely the remnant 
of a binary interacting system, now merged. Following this 
scheme, the cold front may be related to the motion of the 
residual cool core of one of the former interacting clusters, 
with regard to the hot shocked and mixed gas of the present 
cluster. 

In the same way, it is worth noticing that the front in 
Abell 2256 is located around a complex -at least dual- clus- 
ter core, as reveale d by wavelet a nalyses of ROSA T PSPC 
(|Slezak et al.lll994ft and Chandra ([Sun et al.ll2002h images. 
Also in this case, we may infer that the cold front indicates 
the accretion of a subgroup, the core of which is now part 
of the eastern complex system in the cluster. 

8. Neutral hydrogen column density variations 
across the Abell 478 field of view 

The X-ray emission spectra of extragalactic sources are dis- 
torted by the neutral hydrogen absorption along the line of 
sight, whose origin is mostly galactic. We modelled this ef- 
fect by introducing an absorption law which is a function of 
the neutral hydrogen column density, Nh, to the ICM emis- 
sion models n evt (£;, /)F evt (T, Nh, [e]), of Eq. ([2]). Depending 
on the knowledge of the average Nh value in the direction of 
the source observed, and on the spatial variation across the 
source field of view, this parameter may be fixed a priori 
or locally estimated from the observed ICM emission spec- 
tra. The average Nh across the field of view of each cluster 
in our sample was first estimated by fitting Nh absorbed 



emission models to the overall ICM emission spectra. For 
all clusters except Abell 478, these overall spectroscopic 
measurements were consis tent with Nh values measured by 
iDickev fc Lockmanl (|l990l ): for A478, the spectroscopic Nh 
was found to be about twice as high. Due to the inconsis- 
tency seen between both of these values, we left the Nh as 
a free parameter for X-ray spectroscopy. This T — Nh dual 
parameter spectral-fitting process has enabled us to map 
the estimated Nh across the field of view of A478, using a 
similar wavelet algorithm to the ICM temperature mapping 
algorithm described in Sect. 14.1.21 

A map of the spatial distribution of Nh measured by X- 
ray spectroscopy across the field of view of A478 is shown 
by white isocontours on the left panel of Fig. [6j In order 
to reduce a possible spectroscopic bias with brightness gra- 
dient, the map has been obtained from a wavelet analysis 
limited to scales smaller than 2 arcmin. Significant wavelet 
coefficients have been selected according to a hard thresh- 
olding at 2 x a. The Nh distribution appears irregular with 
an east-west elongation, an extended excess to the centre 
of the field of view, and a peak to the to the east. Starting 
from the central excess (Nh ^ 30 x 10 24 m -2 ), the Nh de- 
creases strongly to the north and south of the field of view, 
up to values consis t ent w ith the average galactic value of 
IDickev fc Lockmanl (|l990l ) (Nh < 20 x 10 24 m- 2 ) at a dis- 
tance of 5 arcmin from the centre. 

The Nh central excess sh own on the map of Fig, p 
has already been reported by iPointecouteau et all (|2004t )7 
who performed a spectroscopic analysis within sectors 
of the same data s e t as here. As already proposed by 
IPointecouteau et all (|2004f ), its origin can be investigated 
using the IRAS far-infrared survey of galactic dust emis- 
sion at 100 /im, with angular resolution of about 4 ar- 
cmin. Indeed, the far-infrared galactic dust emission can 
be used as a tracer of the galactic neutral hydrogen column 
density, since dust emission is expected to be correlated 
to the galactic neutral hydrog en at high galac t ic lat itude 
(\b\ > 10 deg), as sho wn by e.g. iBoulanger et al.l (|l996f ) and 
ISchlegel et al.l ((1998). Notice however that this correlation 
is expected to be better constrained within low-Nh regions 
of the sky (Nh < 5 x 10 24 m -2 ) than within high-Nh re- 
gions like the Abell 478 neighbourhood, since the scatter 
of the IR/Nh correlation is higher in such regions, possibly 
du e to the presenc e of m olecular hydrogen, as discussed 
by IBoulanger et all (|l996). A galactic dust emission map 
at 100 jim of the 2 degree neighbourhood of Abell 478 is 
shown in the right panel of Fig. [6j Coming from the last 
gene ration of IRAS maps obtained from the IRIS process- 
ing (Mi ville-Deschenes fc Lagachell2005l ), this map reveals 
the structure of the galactic filament where Abell 478 is lo- 
cated, to a resolution that ca n not b e reached by the radio 
survey of IDickev fc Lockmanl (|l990l ), with angular resolu- 
tion of 1 degree. The IR brightness of this map appears to 
be strongly structured in this sky region, with relative flux 
variations of more than 100 %. Within a 5 arcmin region 
centred on the brightness peak of A478, the IR emissiv- 
ity is of about 18 Mjy/sr. This value corresponds to a Nh 
density column o f 33 x 10 24 m~ 2 fo l lowin g the IR/Nh cor- 
relation factor of IBoulanger et all (p~996). It is consistent 
with density column estimations obtained by X-ray spec- 
troscopy, and encourages us to use the dust emission at 100 
jim as a tracer of the galactic neutral hydrogen distribution 
within this region of the sky. 
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Fig. 6. Left image: IRAS/IRIS 100 jim galactic dust emission map across the field of view of Abell 478 overlaid to the 
neutral hydrogen column density estimated by X-ray spectroscopy (Nh isocontours are equispaced by 2 x 10 24 m -2 and 
decrease from a central value of 30 x 10 24 m -2 ). Right image: IRAS/IRIS 100 fim galactic dust emission map of the 2 
degree neighbourhood of A478, with A478 ICM brightness contours overlaid. The IR emissivity of the black region to 
the south of the cluster is of about 50 Mjy/sr. 



We investigated the accuracy of this tracer by superim- 
posing our Nh map to a portion of the IRAS map in the 
left panel of Fig. [6j Interestingly, we notice some spatial 
correlation between the IR and X-ray estimated Nh map. 
Indeed, both the central extended excess and the eastern 
peak shown on the Nh map correspond to similar structures 
in IR. However, both maps are not fully correlated in flux, 
since the maximum of the Nh map is located to the centre 
of the field of view, while the maximum of the IR map is 
located at the eastern peak. Observations of this sky region 
using infrared data of higher resolution may help to show 
conclusively whether the Nh excess detected near the cen- 
tre of Abell 478 corresponds to a real feature, or if it is due 
to a bias of X-ray spectroscopic estimations. 



9. Discussions and Conclusions 

We present a multi-scale algorithm based on wavelets 
for mapping the ICM temperature structure within clus- 
ters of galax i es. Fo llowing the approach proposed by 
iBourdin et al.l ([2004), the basic scheme of this algorithm is 
i) estimate the searched parameter within square resolution 
elements, sampling the field of view at different scales, ii) fil- 
ter the analysed signal in order to get wavelet coefficients, 
iii) threshold the wavelet transform obtained, in order to 
de-noise the signal and map its 2-D structure. In order to 
perform a more regular reconstruction of structures and to 
reduce thresholding artifacts, we improved the algorithm so 
that, to filter the signal, we use a B 2- spline wave let instead 
of the Haar wavelet used in lBourdin et al.l ((2004) (see para- 
graph [4TL2] for details). This algorithm was used to map 
the X-ray and temperature structure of a nearly complete 
X-ray flux limited cluster sample containing the eight clus- 
ters that, to date, have useful XMM-Newton observations. 
For self-consistency of the results, a similar thresholding 
approach was adopted for each target. 



From previous work, three of the clusters in our sam- 
ple -namely A2029, A1795, and A478- have been identi- 
fied as relaxed systems, while the remaining clusters have 
been identified as major mergers. This classification is sup- 
ported by the power ratios analysis of the surface bright- 
ness of these cluster s, obtained from ROSAT observations 
([Buote fc Tsailll996l ). The XMM-Newton photon images in 
Fig. [1] and the surface brightness contours in Fig. [3] con- 
firm this previous classification: the X-ray morphology of 
the relaxed systems is quite regular and elliptical symmet- 
ric while the morphology of the merging systems is more 
complex with the presence of multiple X-ray peaks. 

In Fig. [3] we show the temperature maps for all clus- 
ters in our sample overlaid to the soft (0.3-2.5 keV) X-ray 
brightness isocontours. Even if with a very different degree 
of complexity, we notice that all clusters, including the most 
relaxed ones, show non radial thermal structure. Let us re- 
call that, since for this work we are only interested in de- 
tecting highly significant thermal structure, we a dopted a 
quite conservative approach which is based on the iDonohol 
(|l995f ) wavelet shrinkage. For this reason the angular res- 
olution of the detected structure appears to be lower than 
the angular resolution expected from EPIC XMM-Newton 
d ata. A wavelet shrinka ge with constant threshold, as tested 
in lBou rdin et al. (2004 ) on simu lated observations, or u sed 
in lBelsole et al.l (|2004l2005l ) and lSauvageot et all (|2005f ) on 
real data, would reveal non radial thermal structure even 
on smaller angular scales. 

Despite the relatively low angular resolution we find 
that, consistent with predictions from numerical hydro N- 
body simulations, the complexity of the thermal structure 
of clusters strongly depends on the dynamical status of clus- 
ters themselves. From Fig. [3] we can see that the temper- 
ature maps of the clusters in our sample may be classified 
into three different types, characterised by an increasing 
regularity of the thermal structure: 
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— "irregular", as for A399 and A401. While selected inde- 
pendently in our sample, these clusters actually form 
a binary and mildly interacting system. It has been 
proposed that the strong irregularity of their tempera- 
ture structure is the result of previous merger activities 
within each cluster of th e system, independently from 
the present interaction ([Sakelliou fc Ponma n 2004). 

— "bimodal temperature", as for A2065 and A2256. These 
clusters are known to show some strong temperature 
anisotropics and to be presently overcoming a late stage 
of merging. They show an elongated geometry of surface 
brightness, with temperature maps characterised by a 
colder and hotter regions on opposite sides along the 
major direction of elongation. This thermal structure is 
dominated by the contrast between cool features to the 
one side, probably associated with accreted material, 
and a hotter region to the other side, possibly assoc iated 
with shocked gas -see e.g. iMarkevitch et all ([1999) and 
ISun et all ((2002), for A2065 and A2256, respectively- 
and separated from the cluster core by a cold front. 

— "regular" , as for the relaxed clusters, A2029, A1795, and 
A478. These clusters have a characteristic cool core re- 
gion surrounded by a hotter gas annulus region peaking 
at about 300-400 kpc from the cluster centre. As with 
the surface brightness, the temperature map of these 
clusters is almost regular and elliptically symmetric, in- 
dicating a good relaxed status for these systems. It is 
worth noting, however, that beyond the overall ellip- 
tical symmetry, we also detect a number of significant 
non radial thermal structures outside their cores. This 
result seems to be consistent with what is predicted by 
hydro N-body simulations of cluster formation, where 
clusters continuously accrete small galaxy groups during 
and between major mergers, which lead to temperature 
irregularities near the region where groups are accreted. 

To conclude, we would like to stress that, despite the 
large dispersion of intrinsic brightness and exposure times 
of the clusters in our sample, the multi-scale algorithm pre- 
sented here enabled us to reveal the thermal structure of 
the ICM of each cluster. Apart the detection of strong fea- 
tures mainly present in merging systems -e.g. cold fronts-, 
we were also able to detect some mild non-radial tempera- 
ture structure outside the core of relaxed clusters. We no- 
tice that the relative temperature variation associated with 
these irregularities at fixed radii is about 10 %. We show 
that these non radial thermal variations affect the measured 
radial temperature profiles (see paragraph Hj) , and lead, at 
least for one of our clusters -A1795-, to an inconsistency 
of mass derivations obtained using hydrostatic equilibrium 
assumption. Nevertheless, for another cluster, A2029, the 
hydrostatic equilibrium assumption appears to be valid de- 
spite the presence of temperature irregularities. 

These results enlighten the possibility of constraining 
the thermalisation status of the ICM and the departure 
from gas hydrostatic equilibrium within bright clusters, 
using data obtained with current X-ray observatories. It 
invites us to further investigations of ICM temperature 
anisotropics within larger samples of relaxed clusters - 
possibly at even larger cluster radii than here by assuming 
specifically an elliptical geometry of the wavelet analysis- 
and to a more accurate evaluation of the dispersions implied 
by these anisotropics on cluster mass estimations. 
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Appendix A: Background modelling for the EPIC 
cameras on board of the XMM-Newton 
satellite 

A typical list of photon impact detection provided by the 
EPIC CCD cameras gathers events associated with various 
processes that have to be taken into account for modelling 
the observed spectra of extended sources. In the case of 
ICM observations, the detected photons are at least related 
to two extended contributions: the observed source itself 
and the cosmic X-ray background (CXB). Furthermore, a 
significant fraction of the events-list is actually associated 
with false detections due to cosmic-ray induced particles 
interacting with the detector. Eventually, a known fraction 
of events is registered during readout periods of detectors, 
which leads to an additive noise due to wrong position reg- 
istration of these so-called "out-of-time" events. 

We model, as an overall "background contribution" to 
the event spectrum n evt (/c,/) F evt (T, Z, Nh, e) of Eq. ([2]), 
a combination of normalised spectral contributions associ- 
ated with CXB, (e), induced particles, F p (e), and readout 
noise, 0(/c, e): 

nbckF bc k(M,e) = E(fc,Z,e) x n C xB (e) 
+n oot (M) 0(fe, Z, e) 

+n p F p (e), (A.l) 

where ncxB, n p and n oot are normalisation terms, 
and where (e) is corrected for spatially variable exposure 
E(fc, Z,e) (see equation [1]), as it is related to a physical 
observation. The CXB spectrum (e) is modelled as the 
combination of a soft radiation associated with foreground 
galactic gas and a broad band contribution accounting for 
extragalacti c bac kgroun d of unresolved AGNs. F ollowing 
iLumb et al.1 (|2QQ2f ) and iKuntz fc Snowdenl ([2000), we use 
a two temperature thermal radiation (fcTl = 0.074 keV, 
kT2 = 0.204 keV) and a power-law (7 = 1.42) for modelling 
the galactic and extragalactic contributions, respectively. 
The cosmic-ray induced particle background F p (e) has been 
modelled from cumulated expositions of telescopes to the 
particles, during in-flight calibration phases. The readout 
noise contribution 0(fc, Z, e) is estimated by integrating and 
normalising the overall signal along CCD columns. While 
normalisation of the readout noise is known a priori for 
a given observing mode, normalisations of both the CXB 
spectrum and particle background must be set by fitting 
our background model, nb c kFb c k(e), to an observed "back- 
ground spectrum". To extract this spectrum, we select 
events from an external ring of the field-of-view with ra- 
dius range of 10.5-12 arcmin, which is a region where back- 
ground emissivity is already dominant but where effective 
area uncertainties are lower than they would be for even 
larger radii. 
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